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' Abstract. We present high resolution {A/AA = 49 000) echelle spectra of the intermediate mass, pre-main sequence stars BF 

^ I Ori, SV Cep, WW Wul and XY Per. The spectra cover the range 3800 ~ 5900 A and monitor the stars on time scales of months 

• 1— ( . and days. All spectra show a large number of Balmer and metallic lines with variable blueshifted and redshifted absorption 

' features superimposed to the photospheric stellar spectra. Synthetic Kurucz models are used to estimate rotational velocities, 

JH I effective temperatures and gravities of the stars. The best photospheric models are subtracted from each observed spectrum 

. . . 1 to determine the variable absorption features due to the circumstellar gas; those features are characterized in terms of their 

velocity, v, dispersion velocity, Av, and residual absorption, iJ^ax- The absorption components detected in each spectrum can be 
grouped by their similar radial velocities and are interpreted as the signature of the dynamical evolution of gaseous clumps with, 
in most cases, solar-like chemical composition. This infalling and outflowing gas has similar properties to the circumstellar gas 
observed in UX Ori, emphasizing the need for detailed theoretical models, probably in the framework of the magnetospheric 
accretion scenario, to understand the complex environment in Herbig Ae (HAe) stars. WW Vul is unusual because, in addition 
to infalling and outflowing gas with properties similar to those observed in the other stars, it shows also transient absorption 
features in metallic lines with no obvious counterparts in the hydrogen lines. This could, in principle, suggest the presence of 
CS gas clouds with enhanced metallicity around WW Vul. The existence of such a metal-rich gas component, however, needs 
to be confirmed by further observations and a more quantitative analysis. 

Key words. Stars: formation - Stars: pre-main sequence - Stars: Circumstellar matter - Lines: profiles - Stars:individual: 
BF Ori, SV Cep, WW Vul, XY Per 
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1. Introduction 
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Observations reveal that the dynamics of the circumstellar (CS) 
gaseous disks around intermediate mass, young main sequence 
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(MS) and pre-main sequence (PMS) stars is extremely com- 
plex. Variable absorption components detected in many lines of 
different elements and ions constitute good examples of such 
complexity. The kinematics and intensity strength of the ab- 
sorption components contain relevant information on the physi- 
cal properties of the gas. Further, their detailed characterization 
and analysis provide clues and constraints on plausible forma- 
tion mechanisms as well as on theoretical scenarios describing 
the structure and nature of the CS gaseous disks. 

The presence of metal-rich planetesimals in the young MS 
P Pic system has been inferred both obs ervationally and the- 
oretically in a long series of papers (e.g. lLagrange et allbOOOt 
and references therein). Summarizing, there are two main argu- 
ments on which this inference is based. Firstly, dust causing the 
far-IR excess and also seen in the (3 Pic disk images may be sec- 
ond generation material continuously replenished by collisions 
of large solid bodies or slow evaporation. Secondly, transient 
spectral line absorptions, usually redshifted, of different chem- 
ical species in a wide range of ionization states can be modelled 
in terms of the evaporation of km-sized bodies on star-grazing 
orbits. Star-grazing planetesimals have also been suggested to 
exist in the 5 1 Oph system, a star with an uncertain evo lution- 
ary status jRoberge et alJbOO^ Ivan den Ancker et alJ l2001): 
also, blueshifted absorption in excited fine structure lines of 
C II* at 1037 A and N ii* at 1085 and 1086 A have been used to 
infer the prese nce of ~1 m bodies in the Vega-type binary MS 
system cr Her dChen & Jura'2003V 

Absorption features similar to those observed in fi Pic have 
been observed towards many HAe stars (see Natta et al. 2000 
for the description of these stars as a subgrou p of HAeBe 
stars) , particularly in the UXOR-subclass (e.g. lOradv et al] 
119961) and, by analogy, they have been interpreted as indica- 
tive of the presence of lar ge solid bodies in the CS disks 
around these PMS stars (e.g. lGradv et al. 2000 . and references 
therein). In principle, this interpretation is not in conflic with 
the accepted time scale for the formation of planetesimals of 
~10'* yr CBeckwith et al.. 2000) . which suggests that planetesi- 
mals should already be present during the PMS phase of stars 
(~1-10 Myr). This explanation for the variable absorption fea- 
tures observedjn HAe s tars is, however, controversial and, in 
fact, ICirinin et alJ ill 994 pointed out other alternatives, more 
concretely dissipation of dust clouds and the simultaneous in- 
fall of cool gas onto the star. Natta et al.. (2000) have ana- 
lyzed the chemical composition of a strong redshifted event in 
UX Ori and shown it to have a solar-like composition. Instead 
of the planetesimal origin for the transient components, those 
authors sugg ested gas accretion from a CS disk. This result is 
supported by Mora et al. ( 2002, from now on Paper I) who pre- 
sented the analysis of a large series of high resolution optical 
spectra of UX Ori. Many variable absorption events in hydro- 
gen and metallic lines were attributed to the dynamical evo- 
lution of gaseous clumps with non metal-rich, roughly solar 
chemical compositions. In addition, iBeust et alJ 1I2OO1I) have 
shown that the p Pic infalling planetesimal model would not 
produce detectable absorptions in typical PMS HAe CS con- 
ditions. We also note that dust disks around HAe stars are pri- 
mordial and can be explained in the context of irradiated PMS 
CS disk models (.Natta et al,200l'l) . 



In this paper we present high resolution spectra of the HAe 
stars BF Ori, SV Cep, WW Vul and XY Per and perform 
an analysis similar to that carried out for UX Ori (Paper it). 
The spectra show very active and complex CS gas in these 
objects; many transient absorption features in hydrogen and 
metallic lines are detected, indicating similar properties of the 
gas around these stars to those of UX Ori CS gas. In addition, 
the spectra of WW Vul show metallic features without obvious 
hydrogen counterparts; in this sense, this star presents a pecu- 
liar behaviour The layout of the paper is as follows: Section|2| 
presents a brief description of the observations. Section |3 
presents the results and an analysis of the photospheric spectra 
and the CS contribution. Section|3presents a short discussion 
on the kinematics and strength of the variable features, and on 
the metallic events detected in WW Vul. Finally, Sect.jsjgives 
some concluding remarks. 

2. Observations 

High resolution echelle spectra were collected with the Utrecht 
Echelle Specti'ograph (UES) at the 4.2 m WHT (La Palma ob- 
servatory) during four observing runs in May, July and October 
1998 and January 1999. 28 spectra were obtained: 4 of BF Ori, 
7 of SV Cep, 8 of WW Vul and 9 of XY Per (XY Per is a 
visual binary sistem with a Herbig Ae primary and a B6Ve 
secondary separated by 1.4", the good seeing during the ob- 
servations allowed us to fully separate both components, in 
this paper only the HAe primary star has been studied). The 
wavelength range was 3800-5900 A and the spectral resolu- 
tion, /l/A/l, = 49000 (6 kms"'). Wavelength calibration was 
performed using Th-Ar arc lamp spectra. Typical errors of the 
wavelength calibration are ~5 times smaller than the spectral 
resolution. The observing log, exposure times and signal to 
noise ratio (SNR) values, measured at /I ^ 4680 A, are given in 
Tabled Further det ails of the observa tions and reduction pro- 
cedure are given bv iMora et alJ 1I2OOII) . For some spectra there 
are simultaneous optical photo-polarimetric and near-IR photo- 
metric observations (Oudmaiier et al. 2001; Eiroa et al ]l200ll) . 
Table[2also presents these simultaneous data. At the time of the 
observations the stars were close to their brightest state, BF Ori 
and XY Per, or at average brightness, SV Ce p and WW Vul 
jHerbst & Shevchenkoll 999tlKiroa et al J2002h. 

3. Analysis of the spectra and results 

3.1. The photospheric spectra 

Circumstellar absorptions with complex profiles and blended 
components are detected in hydrogen and metallic lines in 
all UES spectra of the four stars. The analysis of those spec- 
tral features requires the subtraction of the underlying photo- 
spheric spectra. Such subtrac tion is c arried o ut follow i ng the 
method outlined in detail in IPaper 1 Briefly, iKuruczl (Il993l) 
model atmospheres are used to synthesize photospheric spec- 
tra. Four parameters are estimated: the heliocentric radial ve- 
locity (Vr), the rotation velocity (v sin /), the effective temper- 
atu re (Teff) and the surface gravity (log g). It has been shown 
by iGrinin et aL C2001.) that the photospheric lines of BF Ori 
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Table 1. EXPORT UESAVHT observing logs of BF Ori, 
SV Cep, WW Vul and XY Per. The Julian date (-2450000) 
of each spectrum is given in Col. 1 . Column 2 shows the expo- 
sure time in seconds. Column 3 gives the SNR at A ^ 4680 A. 
Columns 4 to 6 give simultaneous V, %Py and K photopolari- 
metric data, where available. Typical errors are 0.10 in V, 0.05 
in/: and 0.1% in %Pv. 



BF Ori 

Julian date te^p (s) SNR V K 

1112.6324 1800 280 9^65 056 tM 

1113.6515 2700 190 9.79 0.75 7.85 

1209.5542 2700 ITo - ~ 7J6 

1210.4571 2700 280 - 0.14 - 



SV Cep 

Julian date texp (s) SNR V %Pv K 

950.6668 1800 100 - - ~ 

950.6893 1800 120 - 

1025.6260 2700 140 - 0.96 - 

1025.6595 1800 130 - - - 

1026.6684 2700 140 - - - 

1113.4730 2700 70 11.01 1.05 - 

1209.3372 2700 170 - ~ ~ 



WW Vul 



Julian date te^p (s) SNR V %Pv K 



950.6176 


1800 


130 








950.6413 


1800 


130 


10.89 


0.69 




951.6232 


1800 


120 








951.6465 


1800 


150 








1023.5186 


1800 


140 






7.37 


1023.5423 


2700 


190 




0.40 




1112.3689 


1800 


110 


10.77 


0.37 


7.44 


1113.3958 


2700 


120 


11.03 


0.65 


7.50 




XYPer 


Julian date 




SNR 


V 


%Pv 


K 


1024.6728 


1800 


200 








1024.6967 


1800 


170 


9.04 


1.49 




1025.6948 


1800 


230 




1.55 




1025.7171 


1800 


270 








1026.7065 


1800 


220 








1112.4978 


1800 


270 


9.12 


1.65 


5.97 


1113.5483 


2700 


230 


9.05 


1.53 


5.99 


1207.3786 


2700 


190 








1209.3844 


1800 


260 


9.51 


1.58 


6.18 



and WW Vul can be well reproduced using solar metallicity 
synthetic spectra. On the other hand, the presence of CS com- 
ponents in most of the lines makes it very difficult to carry 
out a detailed abundance analysis. We have thus decided to as- 
sume that the four stars have solar metallicities. Th e atomic line 
data h ave been obtained from the VALD database ( iKupka et alJ 
Il999h . The 'best' synthetic model, defined by the parameters 
Usted before, is selected by comparing some appropriate faint 
photospheric absorption lines among the observed spectra and 
the synthetic ones. This is not straightforward because each 
star behaves differently, and the choice of pure photospheric 



Table 2. Stellar parameters defining the 'best' synthetic Kurucz 
models for each star. See Sect. 13. II for a discussion on the un- 
certainties of TefF, log g and v sin /. 



Star 


ieff 

(K) 


logg 


V sin I 
(km s-') 


v,. 

(kms-') 


BF Ori 


8750 


3.5 


37 


23.1 + 1.9 


SV Cep 


10000 


4.0 


225 


-11.9 ±0.8 


WW Vul 


9000 


4.0 


210 


-104 ± 1.2 


XYPer 


8500 


3.5 


200 


8.3 ± 0.6 



lines in the spectra of such highly variable objects is not triv- 
ial. Therefore, slightly different, ad-hoc approaches for each 
star are needed. These approaches are discussed below and the 
stellar parameters giving the best synthetic spectra are given in 
Tabled They are compatible within the uncertaint ies with the 
spectra l types and rotational velocities quoted b y iMora et all 
(2001) and with the results bv lOrinin et alJ ( l200ll) . who studied 
BF Ori and WW Vul. 

BF Ori: The photospheric lines are narrow and do not show 
a noticeable variability, v^ is estimated using faint photospheric 
lines and its rms error is low. Estimated errors are 250 K for Teif 
(the step in the synthetic spectrum grid created using Kuruzc's 
models), and 6 kms"' for vsin/ (the spectral resolution); the 
log g values considered have been restricted to 3.5 and 4.0. 
Fig. [2 shows the excellent agreement between the synthetic 
model and the observed median spectrum. The extra absorption 
seen in the stronger lines is due to the circumstellar contribu- 
tion. 

SVCep: According to lFinkenzeller & .Tankov ics'('l984'') the 
radial velocity of HAe stars coincides within a few km s ' with 
the radial velocity of interstellar (IS) lines. The strong Nai D 
IS components in the line of sight of SV Cep allow us to de- 
termine a precise value of the radial velocity of this star Its 
high temperature, rotation velocity and variability make it ex- 
tremely difficult to identify photospheric lines in order to de- 
rive the Kurucz stellar parameters. We have circumvented this 
pro blem by using the EXPORT intermediate resolution spec- 
tra (" Mora et al]l200l V which have very high quality Ho- pro- 
files and are of great help in making a first estimate of Tetr 
and log g. ^From the Ha wings and the almost complete ab- 
sence of weak photospheric lines in the high resolution spec- 
tra, Teff - 10000 K and log g = 4.0 are estimate d. It was 
assumed that vsin/ = 206 kms-' (lMoraetalJl200ll) . The ad- 
equacy of this choice can be seen in Fig. |2] where the ob- 
served median SV Cep Ha line and 2 synthetic Balmer pro- 
files (Teff = 10000 K, log g = 4.0, vsin / = 206 kms"' and 
Teff = 10000 K, log g - 3.5, vsin / = 206 kms"') are shown. 
Using those values we have identified 7 blends of weak pho- 
tospheric lines (absorption intensities < 4% of the contin- 
uum) with very low variability. The blends are at ~ 3913, 
4129, 4176, 4314, 4534, 5041 and 5056 A. A grid of pho- 
tospheric spectra has been generated with different values of 
Teff and vsin/ (log g was assumed to be 4.0 from the Ha 
analysis) and the rms differences between the synthetic and 
observed spectra have been estimated. The lowest rms differ- 
ence is obtained for Teff - 10000 K and vsin / = 225 kms"'. 
The differences between the Ha wing profiles broadened to 
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Fig. 1. Synthetic (dashed line) vs observed spectrum (solid line) for BF Ori. The agreement between the photospheric and the 
observed median spectrum is remarkable over a large wavelength range. The extra absorption in the stronger lines (identified in 
the figure) is due to the circumstellar contribution. This example illustrates the need to use very faint lines to estimate the stellar 
parameters. (This figure is available in color in electronic form) 



vsin/ = 225 kms"' and vsin/ = 206 kms"' are negligible, so 
it was not needed to compute log g again. Fig.|2lalso shows the 
best synthetic and the observed median spectra for the selected 
blends. Uncertainties of -500 K (the step of the Kurucz's mod- 
els) and ~IQ% for Teff and v sin / , respectively, are estimated. 

WW Vul: V,- is estimated from the sharp Na i D IS lines. The 
star was very active during the observing runs and its spectra 
show a large number of broad and variable absorption features, 
e.g. redshifted features are superimposed on practically each 
photospheric line in the JD 1113.40 spectrum, which pose se- 
vere difficulties for the selection of appropriate photospheric 
lines. Nevertheless vsini has been estimated with the Mgii 
4481 A blended doublet. Many pairs of Teu - log g values can 
reproduce the observed spectra, though the agreement between 
the observed and synthetic spectra is generally poor. We have 
selected Tgff = 9000 K, log g - 4.0 as representative of the 
WW Vul photosphere because of its compatibility with the re- 
sults by Moraet al. (2001), though large errors are likely, but 
not larger than around 10% in both Teff and v sin /. Fig.|3]shows 
a comparison between the observed median spectrum and the 
broadened synthetic one in two spectral regions where a good 
fit is achieved; this figure also shows the unbroadened Kurucz's 
model with some line identifications. 

XY Per: v,. is estimated from the Nai D IS components. 
The photospheric lines are very broad but the high SNR of the 
spectra and the relatively low Teff allow us to identify a large 
number of faint line blends in order to perform a precise esti- 
mate of Teff, log g and vsin/ (16 faint blends with absorption 
intensities < 8% of the continuum with very little CS activity 
could be identified). Errors of Teff and log g are of the order 




■ Observed spectrum — 

■ Synthetic spectrum 

■ Synthetic spectrum, vsini - 



4200 4210 4220 

Wav elength (A) 



Fig. 3. The observed median UES spectrum of WW Vul (black 
continuous line) compared to the broadened synthetic one (blue 
dashed lines) in two different spectral regions. At the top of the 
figure, the unbroadened Kurucz model is shown with some line 
identifications (red solid line). (This figure is available in color 
in electronic form) 



of the step in the Kurucz models, 250 K and 0.5 respectively, 
while the error in vsin / is very low, <10%. The comparison 
between the synthetic and the observed median spectra of the 
16 blends is shown in Fig.|3 
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Fig. 2. Synthetic vs observed spectra for S V Cep. Left: comparison of the intermediate resolution Ho' spectrum and two syn- 
thetic Wa profiles. The b est result is obtained for Teff - 10000 K and log g - 4.0 (blue dashed line). It was assumed that 
vsm; = 206 kms-i jMoraet al.ll200ll) . Right: Median UES spectra (solid lines) of several blends compared with the best syn- 
thetic spectra, with Tgff - 10 000 K, v sin / = 225 km s"' (blue dashed lines). (This figure is available in color in electronic form) 
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Fig. 4. Synthetic (dashed lines) vs observed median (continuous lines) spectra of XY Per for 16 spectral features with very 
low CS contribution. The stellar parameters of the synthetic spectrum are given in Table 121 (This figure is available in color in 
electronic form) 



3.2. The circumstellar transient absorption contribution 

Once the best photospheric spectrum of each star is determined, 
the circumstellar contribution to each observed spectrum can be 
estimated by subtracting the synthetic one. The residual spec- 
tra show transient absorption features, which can be character- 
ized by means of t he normalized resi dual absorption, defined as 
R^l- Fohs/Fsyn jNatta et al.l2000l) . The profile of each line 
reflects the blending of several components. A multigaussian fit 
providing the radial velocity, velocity dispersion and absorp- 
tion stre ngth is used to identify the individual components (see 
IPaper iL for details). Broad redshifted and blueshifted absorp- 
tions at different radial velocities are found in the Balmer and 
metallic lines for all 4 stars analyzed. We apply the multigaus- 
sian fit to Balmer lines (H/3 4861 A, Hy 4340 A, 4102 A, He 
3970 A, 3889 A), Can K 3934 A, Can H 3968 A, Nai D2 



5890 A and Nai Dl 5896 A, as well as to fainter metallic 
lines Fell 42 multiplet (a6S-z6Po triplet: 4924 A, 5018 A and 
5 169 A), Ti 11 4444 A, Ti ii 4572 A, Fe i 4046 A, Sc ii 4247 A 
and Ca 1 4227 A. We have chosen these ionic lines because they 
show significant CS variability and are relatively strong and 
isolated. Narrow IS components (mainly Na i, Ca ii and Fe ii) 
with the stellar radial velocity are also detected for all the stars. 

Consecutive spectra with a time delay of ~1 hour of 
SV Cep, XY Per and WW Vul were taken on several nights 
(see Table 1). These spectra were quite similar and the gaussian 
deconvolution of the R profiles essentially provides the same 
values for the fit parameters; thus, any significant variation of 
the phenomena causing the transient absorptions is excluded 
on this time scale, at least during these observing periods. This 
result gives us confidence in the identification of the compo- 
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nents and allows us to add the spectra taken during the same 
night in order to increase the SNR. Tables |6l0|8] and |9] avail- 
able at the CDS, give the radial velocity shift v, the velocity 
dispersion Av and the absorption strength Rmax, the peak of the 
R profile, of each identified broad transient absorption compo- 
nent of the lines listed above for BF Ori, SV Cep, WW Vul and 
XY Per, respectively. Column 1 gives the the corresponding 
Balmer or metallic line. Col. 2 gives the Julian Date, Col. 3 rep- 
resents the event assigned to the particular absorptions. Col. 4 
gives the radial velocity shift v. Col. 5 lists the velocity dis- 
persion Av and Col. 6 gives the absorption strength Rm-^x- JD 
values in Tables and |5] correspond to the starting time of 
the first spectrum of each night. In the following, whenever a 
JD is given, 2450000 is subtracted. 

Tables |6j |S] and |5] show that absorption components 
with similar radial velocities appear/dis appear simultaneously 
in different lines, as is observed in UX Ori (iPape^. We as- 
sume that absorptions with similar velocities come from the 
same gaseous clump, which can be characterized by an aver- 
age radial velocity <v> (a Transient Absorption Component or 
TAC). The time evolution of the TACs' velocity is referred to as 
an event and represents the dynamical evolution of the gaseous 
clumps. We point out that there is an uncertainty in identify- 
ing TACs detected on different nights with the same g aseous 
clump; our assumption relies on the UX Ori results (.Paper iL 
see also below). Figs. Bl 1711^ and fTTI plot <v> of the identified 
TAC versus JD. <v> is a weighted average in which the lines 
H/3, Hy, U6, Ca ii K, Fe ii 4924 A and Fe ii 50 1 8 A have an ar- 
bitrarily assigned weight of 1 because of their higher intensity 
and non-blended nature. A weight of 1/2 is also arbitrarily as- 
signed to the remaining lines in order to reflect that they are ei- 
ther blended (He, Ca ii H and Fe ii 5 169 A), weak (Sc ii 4247 A, 
Fei 4046A and Cai 4227 A) or affected by telluric lines (Nai 
D2 & Dl). The weighted number of lines used in each average, 
which can be fractional because of the 1/2 weights, is plotted 
next to each point. The TACs have been grouped according to 
the event they represent; thus, the figures show the dynamical 
evolution of the gaseous clumps. Note that some events were 
only detected once (only one TAC). Figs.l6ll8lll0land[T2lshow 
the R profiles of some selected lines for the four stars. The 
line absorption components and the corresponding event iden- 
tifications are indicated. The lines are H/3, Can K (except for 
BF Ori), Nai D2, Nai Dl and Feii 5018 A. Hy is shown for 
WW Vul (Fig. llOI I since H/3 has a large underlying emission. 
The results for each star are presented in the following. 

BF Ori: 5 TACs grouped in 3 different events were detected 
in the spectra of BF Ori (Fig. |5}- Event #1 is an accelerating 
redshifted event, first detected in the JD 1112.63 spectrum at 
approximately the stellar radial velocity. Absorbing gas is seen 
in Balmer and metallic lines of Na i, Fe n, Ti ii, Sc ii, and Ca i. 
The strongest H i lines (Hy6, Hy and HS) appear to be saturated, 
i.e. Rmax is very close to unity, and the metallic lines are also 
very strong. The Balmer lines are broader than the metallic 
ones. The parameters of at least some of these lines might be 
influenced by IS gas absorption (a careful look at the Nai D 
lines shows the presence of two peaks in the JD 1 1 13.65 spec- 
trum. Fig. |5Jl. The velocity dispersion of the lines tends to be 
larger when the event increases its velocity, while R^ax val- 



ues tend to decrease, though changes are modest. Event #2 is a 
strong redshifted decelerating event; the behaviour of its TACs 
is in general similar to those of #1. #3 represents blueshifted 
gas only detected in the last night of January 99 and is fainter 
than the redshifted ones. The R^ax values of the Balmer lines 
are low, but they might be saturated, because Rmax does not ap- 
pear to decrease as we follow the series in what will be called 
the "expected Balmer decreasing trend" for optically thin gas 
(R^,x(m > Rmax(iiy) > R^,x(m > /?max(HO); also the rel- 
ative intensity of the metallic lines with respect to the H i ones 
seems to be lower 

The Can K line has non-photospheric profiles with the 
simultaneous presence of redshifted and blueshifted compo- 
nents, but their radial velocities do not match the absorp- 
tions observed in other lines (except the blueshifted TAC in 
JD 1210.45). 

SV Cep: Broad absorptions of H i, Ca ii, Na i and Fe ii are 
detected in the spectra of SV Cep, but no variability is found 
in Ti II, Sc ii, Fe i and Cai (unlike the other stars in the paper). 
The broad absorptions represent 10 TACs grouped in 8 different 
events: 3 of them correspond to outflowing gas and the remain- 
ing 5 to infalling gas (Figs. Q and Blueshifted gas shows 
small radial velocities, on average < 20 kms"', while red- 
shifted components display velocities as high as 160 kms"'. 
Only one spectrum was taken in May 98, October 98 and 
January 99, i.e. the events of these periods are composed of 1 
TAC only: these data represent isolated snapshots of the CS gas 
around SV Cep and no temporal evolution can be inferred from 
them. The 4 TACs detected in July 98 can be grouped in two 
events: #3 corresponds to redshifted gas with practically con- 
stant radial velocity, and #4 is gas observed at a velocity close 
to the stellar radial velocity or slightly blueshifted. In general, 
the /?max values of the Balmer lines show the expected Balmer 
decreasing trend and are broader and much stronger than the 
metallic lines, ^max(H5)/^max(Fe ii) > 3. There are, however, 
some exceptions. #1 shows relatively strong Fe ii and it is not 
clear that the Balmer lines show the expected decreasing trend. 
In #5 the H i lines are not broader than the Fe ii ones, and the 
Balmer lines are probably saturated (but note that this event is 
very faint, and it could be an overinterpretation of the fit proce- 
dure). The strongest H i lines are saturated in #7 and also in the 
JD 1026.66 spectrum of #3. There are anticorrelated changes 
in the Av and Rm^x values of #3, but they show the same trend 
in #4. 

WW Vul: 15 TACs grouped in 9 different events are iden- 
tified in the 5 spectra of WW Vul (Figs. l9l and I10>. 4 events 
are seen in May 98. #1 is redshifted gas with saturated Balmer 
lines and strong Fe ii lines. Both H i and Fe ii lines have similar 
Av and from one TAC to another R^^x and Av show opposite 
trends. #2 corresponds to low velocity blueshifted gas clearly 
detected in the metallic lines but no counterpart in the hydro- 
gen lines is apparent (see the Feii 5018 A line in Fig. llOt . #3 
and #4 are blueshifted accelerating events. The only spectrum 
of July 98 (JD 1023.51) reveals 3 TACs: #5 is a broad, Av > 
100 km s"' , redshifted event only detected in metallic lines (the 
broad wing of the Feii 5018 A line profile in Fig.llOL #6 is a 
very low velocity blueshifted component (v ^ -5 kms"'); this 
event is significantly broader in the Hi fines (which are sat- 
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Fig. 5. Events in BF Ori. Each point coiTesponds to the radial velocity of one TAC and represents the average velocity <v> of 
the absorptions with similar radial velocities detected in different lines. <v> is a weighted average in which the higher intensity 
non-blended lines are assigned a weight of 1 and the rest 1/2. Error bars show the rms error of the average velocity; the numbers 
above the data points indicate the weighted number of lines used to estimate the average. Fractional numbers arise from the 1/2 
weight attributed to some lines (see text). Redshifted events (infalling gas) are printed in red colour, while blueshifted events 
(outflowing gas) are in blue. (This figure is available in color in electronic form) 



urared) than in the metallic ones (the IS contribution cannot 
be separated from this low velocity event). #7 is a relatively 
narrow blueshifted event only detected in metallic lines and 
is clearly distinguished as a peak in the line profiles. Again, 
metallic redshifted absorptions at v^90kms ' without Hi 
counterpart are detected on JD 1112.36 (October 98). Similar 
metallic absorptions are also detected on JD 1113.39, but on 
this date saturated hydrogen components with basically the 
same kinematic parameters (including the velocity dispersion) 
are present. We tentatively identify both TACs with the same 
event, #8, although we cannot exclude the possibility that the 
metallic absorptions detected on each date could be due to dif- 
ferent gas. Finally, #9 is a strong, redshifted, decelerating event, 
identified in H i and in many metallic lines. In this case, the hy- 
drogen lines are considerably broader than the metallic ones. 
The Hi lines seem saturated on the first night, while on the 
second one the expected Balmer decreasing trend is observable 
and on both nights the Fe ii lines are relatively strong. 

XY Per: 16 TACs grouped in 9 different events are detected 
in the 7 XY Per spectra. #1 is a faint redshifted event detected 
on the last two (out of three) July 98 nights. #2 corresponds 
to blueshifted accelerating gas; all three TACs of this event are 
strong and broad. #3 corresponds to a relatively narrow and 
faint blueshifted event only detected on JD 1025.69. #4 is also a 
narrow and faint event, centered at approximately the stellar ra- 
dial velocity, only detected on JD 1026.71; IS absorption could 
contribute to the metallic Rmax values. #5 represents decelerat- 
ing redshifted gas. #6 is a strong redshifted event only detected 



on JD 1 1 12.49. #7 corresponds to decelerating blueshifted gas; 
the expected Balmer decreasing trend is clearly present in both 
TACs and metallic lines are very strong. Both Av and Rmia in- 
crease very significantly from JD 1112.49 to JD 1113.55. #8 
is decelerating redshifted gas detected on January 98 (note that 
the time interval between the two TACS grouped in this event 
is 48 hours); in both TACs the metallic lines are very broad and 
strong: Fe ii lines are even broader and stronger than the Balmer 
ones which are faint and seem saturated. Finally, #9 represents 
accelerating blueshifted gas. 

4. Discussion 

The objects studied in this work and UX Ori are very 
much alike. All are PMS HAe stars. UX Ori ^ BF Ori an d 
WW Vul are bonafide UXOR-type objects (e.g. Grinin'2000|), 
i.e. their light curves show high-amplitude variability (Am > 
2.0 mag), Algol-like minima, a blueing efifect and an in- 
crease of the polarization when the object brightness d ecreases. 
SV Cep also shows UXOR characteristics (Rostopchina et alJ 
bOOOtlOudmaiier et alJEool . Further data are required before 
XY Per can be confidently classed as a UXOR-type object 
( Oudmaiier et al. 200D . although it does share the same com- 
plex and variable spectroscopic behaviour of the other stars. 

The time covered by the UX Ori spectra allowed us to an- 
alyze and identify the TACs as due to the dynamical evolution 
of gaseous moving clumps; this identification and dynamical 
evolution was particularly convincing in the case of a clump 
detected in four spectra taken within a time interval of four 
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Fig. 6. BF Ori R profiles. The normalized residual absorption profiles (/?=!- Fobs / Fsyn) of Hj3 (left), Na i D2 and Dl (middle) 
and Fen 5018 A (right) are shown in the figure (black colour). The corresponding spectra are indicated in the vertical right axis 
(Julian and civilian epochs). The identified TACs (gaussian components, pink colour) and the reconstructed R profile fit (blue 
colour) are displayed. Event numbers are shown under the gaussians. The zero velocity interstellar components are marked as 
"IS". (This figure is available in color in electronic form) 
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hours fraper f). We rely on those results to ascribe to the same 
gaseous clump groups of H i and metallic absorption compo- 
nents described in the previous Sect. We are aware, however, 
that this identification is doubtful in some cases, since the num- 
ber of spectra for the stars in this paper is smaller and the time 
coverage is poorer than for UX Ori. With this cav eat in mi nd, 
we will discuss these spectra as we have done in IPaper ll for 
UX Ori. We are quite confident that the main conclusions of 
this paper are not affected by the uncertainties with which some 
specific event can be identified. 

4.1. Kinematics 

Accelerating/decelerating blueshifted and redshifted events 
were detected in UX Ori, the events seemed to last for a few 
days and their acceleration rates were a fraction of a ms"^. 
Within the present time coverage limitations, the same is ob- 
served in the outflowing and infalling gas of BF Ori, WW Vul, 
SV Cep and XY Per. All 5 stars share the trend that infalling 
gas shows the largest velocities (the exception is event #4 of 
WW Vul in the JD 951.62 spectrum) and that blueshifted ab- 
sorptions are detected when redshifted absorptions are present 
(the exception might be the JD 1024.67 spectrum of XY Per 
which shows a low velocity outflow but no in falling gas). 
Simila r results are also present in the spectra of lOrinin et alJ 
il200ll) . 

Infalling gas appears to have larger dispersion velocities 
than outflowing gas. In addition, the Hi lines are broader than 
the metallic ones in approximately 40 % of the identified TACs, 



as indicated, for example, by the fact that Avpen < 0.66Avhi (in 
UX Ori this trend was noticed in the Nai lines). The fraction 
of TACs with broader H i lines varies from star to star. For ex- 
ample, in 12 out of 16 detected TACs in XY Per Hi and Feii 
have similar velocity dispersion, and only 3 out of 10 TACs in 
S V Sep follow this trend. On the other hand, the Fe ii lines are 
broader than the Hi lines in only two events, #5 of SV Cep 
and #8 of XY Per. It is worth noting that the event of XY Per 
has metallic /?max values at least as large as those of the Balmer 
lines and that the S V Cep event is very weak (see above for this 
event). 

There seems to be a correlation between the dispersion ve- 
locity and the velocity of the TACs in the sense that TACs ap- 
pear to be broader when the velocity increases. If we perform 
a linear regression (Av = A + B x |v|) to the whole set of data, 
we find a correlation coefficient of 0.68. There is also a sugges- 
tion of an anticorrelation between <Av'Fen> and the < Rm^x > 
values of the TACs during their evolution, i.e. events become 
fainter when they increase their dispersion velocity. This was 
also suggested in the case of UX Ori. For every event with 
more than 1 observation (14 in total), <AvFen> and < Rma^ > 
have been normalized to the first observed values. If 2 'anoma- 
lous' events, which are the most uncertain identifications in the 
whole sample, are removed (WW Vul #8 - one TAC is ob- 
served only in the metallic lines while the second TAC is de- 
tected in both metallic and H i lines -, and XY Per #7 - the 
Av and < Rmax > variations are much more extreme than those 
of all other detected events) we find a linear correlation coeffi- 
cient of -0.66. Both correlations are highly significant, i.e., the 
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Fig. 9. Events in WW Vul. Events marked with a "#" ( 1 , 3, 4, 6, 8 and 9; red for redshifted events and blue for blueshifted events) 
are detected both in hydrogen and metallic lines. Events denoted with an "@" (2, 5 and 7; green) are only seen in metallic lines. 
The square point of #8 in JD 11 12.37 correspond to TACs only observed in Fe ii and Ti ii. Further details as for Fig|5] (This figure 
is available in color in electronic form). 



probability of randomly obtaining such coefficients from 2 un- 
related variables is <1% for the data sets considered. We have 
to point out, however, that the <Av'Fen> vs < Rmax > anticorre- 
lation dissapears if the two 'anomalous' events are included in 
the statistics. 

4.2. Line intensity ratios 

Many absorption components in Tables |6l to O have R^ax val- 
ues close to unity, which suggests that they are saturated. 
This is not the case even for the strongest events of UX Ori 
(Paper 1) (in fact, the ratio among the Rm^^ line values of the 
24 UX Ori TACs does not vary much, which allowed us to 
estimate a line residual absorption average). Nevertheless, we 
have followed the procedure of Paper 1 to investigate whether 
a "fixed" ratio among the line absorption strengths might be 
present in those TACs which are most likely unsaturated. Thus, 
we have excluded from this exercise lines wich show signs 
of being saturated, i.e. those with Rmax > 0.8 and TACs with 
Balmer lines of similar strength. The line Feii 4924 A has 
been taken as a reference as it shows the lowest statisti- 
cal errors (other lines which have been considered are H^, 
Can K and Feii 5018 A). The ratio max, line /^max.ref. line (e.g. 
^max,H5/^max,Feii4924A) ^as been Computed for each line of ev- 
ery TAC, and later the mean </?max,iine/^max,ref.iine> has been 
estimated using a sigma-clipping algorithm to reject bad points. 
Table 13 gives </?max, line /^max.ref. line > for cach Star, together 
with statistical errors and line rejections (%). 

Values quoted in Table |3l can be used to assess whether 
the gas causing the variable absorptions is optically thin or 



not, simply by comparing the Rma^ ratios of lines belonging 
to the same element multiplet with their gf ratios (we recall 
that the TACs suspected of being saturated have already been 
excluded). We have analysed the Balmer, the Nai D doublet 
and the Fe ii 42 Unes, for which H6, Nai D2 and Fe n 5018 A 
have been taken as reference, respectively. Rmax ratios have 
been computed following the above procedure and gf values 
have been taken from the VALD database ( Kupka et all l 19991) 
for Hi and Nai. Feii lines do not have reliable experimen- 
tal gf values, partly because Feii 5169 A is blended with 
Mg I 5167 A, Fe I 5167 A and Mg I 5173 A. Following the 
suggestion made by T.A. Ryabchikova and E.G. Kupka (private 
communication), we have used the semiempirical values com- 
puted according to Raassen & Uvlings ( 1998) and available 
at http://www. science.uva.nl/research/atom/levels/levtext. html 
The comparison between the estimated ratios and the theoret- 
ical gf shows that the Fe ii 42 triplet is most likely saturated 
in BF Ori but not in the rest of the stars, and that the Nai D 
doublet is also probably saturated in BF Ori and WW Vul 
(Table 0}. Concerning the Balmer lines (note we are referring 
to those TACs apparently unsaturated) their ratios are very 
different from the gf ratios. This is, in principle, similar to the 
case of UX Ori wh ere the lines do not seem to be saturated 
in any of its events. Paper? suggested that the UX Ori results 
could be explained by underlying line emission caused by a 
spherical occulting cloud with a temperature Tex ~ 7000 K 
and a radius of the order of the UX Ori corotation radius, 
Rc\oud/R, ~ 1-6 (see iRodgers et alJl2002l for details of the 
assumptions). However, only XY Per presents Balmer line 
ratios which could be adjusted using this scenario, namely 
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Fig. 11. Events in XY Per. Details as for Figls] (This figure is available in color in electronic form). 



gas at a Tex ~ 6600 K at approximately the corotation radius 4.3. Origin of the variable circumstellar gas clumps 
Rcioud/R* ~ 1-6. This is not appHcable for the rest of the detected in hydrogen and metallic lines 

stars, where it is most likely that the Balmer lines are always 
saturated. 

Excluding the fact that redshifted events seem to be detected at 
larger velocities than blueshifted events (a fact deserving fur- 
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Table 3. Ratios of the average R^nx parameter of several lines to Feii 4924 A < /^max.iine/^max Fen4924A ^) ^'^^ ssich star. The 
values correspond to lines in TACs which are most likely not saturated (see text). Values with no error mean that only one TAC 
is available. The sigma-clipping threshold adopted is 2.0 cr, except for the values followed by the symbol t, in which 1.5 cr has 
been used. The percentage of rejected lines is given in brackets. 
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5.23 


± 1.82 
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±0.83(14%) 
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± 2.29 


4.53 
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±0.82(11%) 


HS 






4.96 


± 1.85 


4.83 


±0.97 


2.86 


± 1.15 


Hf 


1.50 ± 


0.30 


3.71 


± 1.27 


2.17 


±0.69 


2.09 


±0.66 


CanK 


5.79 




4.05 


±0.97(10%) 


2.37 


±0.65 


2.20 


± 1.16(8%) 


NaiD2 


1.06 ± 


0.40 


1.27 


±0.40 


1.45 


±0.64(17%)'" 


0.76 


±0.12(20%) 


NaiDl 


0.89 ± 


0.39 


0.85 


±0.33 


1.50 


±0.92 


0.48 


± 0.09 (11%) 


Fen5018A 


1.08 ± 


0.15 


1.41 


± 0.14 


1.28 


±0.11 (15%) 


1.16 


± 0.17 (14%) 


Fen5169A 


1.18 ± 


0.15 
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±0.30 
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±0.13 (15%) 
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± 0.25 (7%) 
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0.57 ± 


0.29 
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0.45 ± 
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0.47 ± 
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Ca I 4227 A 


0.32 ± 


0.06 






0.18 
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Table 4. Estimated < /?max, line /^max.ref. line > I'atios among the lines of the Nai D doublet and the Fe ii 42 triplet for each star The 
theoretical ratios (gfune / gfref.iine) and the line taken as reference are given in the last Col. The Fe ii 5169 A line is included for 
comparison purposes though its ratios are likely affected by the blend with Mg 1 5 167 A, Fe i 5 167 A and Mg i 5 173 A. 



Line 


BF Ori 


SVCep 


WW Vul 


XY Per 


Theor. 


Reference 


NaiDl 


0.83+0.11 


0.73+0.15 


0.81 ±0.12 


0.63 ± 0.07 


0.50 


NaiD2 


Feii4924A 


0.94 + 0.16 


0.72 ± 0.08 


0.79 ± 0.07 


0.86 ±0.17 


0.69 


Fen5018A 


Pen 5 169 A 


1.09 ±0.03 


1.07 ±0.15 


1.06 ±0.13 


1.07 ±0.19 


1.25 


Fen5018A 



ther observations and also a theoretical explanation), the gas is 
observed in the same absorption lines with similar R and veloc- 
ity dipsersions, i.e., there are no fundamental differences in the 
behaviour of the infalling and outflowing gas in all the observed 
stars, including UX Ori. This result suggests that their physical 
conditions are rather similar and that they probably originate at 
similar distances from the star Roughly solar abundances were 
found by .Natta et al. (2000.) for an event of redsh ifted gas in 
UX Ori, which was further supported by IPaoer ll This led to 
the conclusion that the clumps of CS gas in UX Ori are non 
metal-rich. Since high velocity gas is observed simultaneously 
in the Balmer and metallic lines in BF Ori, SV Sep, XY Per 
and in most of the detected events in WW Vul, a similar con- 
clusion very likely holds also for the CS gas in these stars (but 
see next Sect, for the metallic events in WW Vul). Unlike the 
detected events in UX Ori, the gas in those stars is often op- 
tically thick, as suggested by the saturation of the absorption 
features. This might indicate that high density gas is more fre- 
quently observed in BF Ori, S V Sep, WW Vul and XY Per than 
in UX Ori. 

IPaper j compares the dynamics of the gaseous clumps in 
UX Ori with the predictions of magnetospheric accretion mod- 
els (see Hartmann 199Si for a very good ba sic d escription 
of thi s theory) and diff erent wind models (e.g. lOoo dson et al. 
Il997t Ishu et aUKOOd iKonigl & Pudrit^bOOff ). The present 
data do not add any new substantial aspect to that discussion, 
only that, with very few exceptions, outflowing gas displays 
smaller velocities than infalling gas. Thus, to avoid repetition 
we refer to that paper, stressing the need for further theoreti- 
cal efforts to explain the complex circumstellar environment of 
HAe stars, at least to the level of understanding achieved for 
the less massive T Tauri stars. 

4.4. The intriguing case of WW Vul 

WW Vul is a very interesting case, somewhat different from 
the other stars. In all the spectra we have obtained, in addi- 
tion to events detected in both metallic and hydrogen lines, as 
it is always the case in the other objects, we see also metal- 
lic absorption features (both blueshifted and redshifted) that do 
not seem to have a counterpart in the hydrogen lines. They 
appear as broad high velocity wings in the R profiles, e.g. 
events #5 and #8 in the JD 1023.52 and 1112.37 spectra, as 
well as relatively narrow distinct peaks, e.g. #7 in JD 1023.52 
(Fig. not . The observed outflowing events are narrower and 
have larger /^max values than the infalling ones. To ensure that 
the metallic events are not artificially introduced by the multi- 
gaussian fit we have proceeded in two different ways: 1 . We 



have carefully inspected the original spectra, i.e., before sub- 
traction of the photospheric components. Fig.^jshows several 
lines as observed in the original (without photospheric subtrac- 
tion) JD 1023.52 WW Vul spectrum. In particular, the relatively 
narrow, blueshifted, metallic event #7 is clearly seen on a sim- 
ple visual inspection as an absorption peak with no distinctly 
separate hydrogen counterpart. We do note that there is hydro- 
gen absorption covering the velocity range of event #7 and that 
the hydrogen lines are quite broad and likely to be saturated at 
this time. As a result, Gaussian fitting may not properly repre- 
sent the intrinsic shape of the hydrogen components; however, 
without advance knowlegde of their intrinsic shape, Gaussian 
fitting is used here as a default procedure. 2. In general, our ap- 
proach has been to use a number as small as possible of gaus- 
sians to reproduce the R profiles in a self-consistent way, avoid- 
ing to overfit the data. 

Nevertheless, since we were intrigued by the WW Vul be- 
haviour, we have performed a number of additional numeri- 
cal tests, to check as well as possible that the apparent lack 
of a hydrogen counterpart to some event seen in the metallic 
lines was not an artifact. First of all, we have added to the H i 
identified components more gaussians which take into account 
the parameters of the absorptions only identifed in the metallic 
lines. This led to unphysical results, i.e., adding more gaus- 
sians to the H I lines does not produce kinematic components 
similar to those observed in the metallic lines, but create spu- 
rious Hi components. For example, in JD 1023.52 we tried to 
fit 3 components to the Hi /?niax profiles, i.e. the same num- 
ber of gaussians obtained from the Fe n fit. The data obtained 
from the Balmer lines deconvolution are not self-consistent: i) 
strange, unrealistic Balmer intensity sequences are obtained, 
e.g. in event #7 Rn,^x.H( = ^max.H/j > ^max.Hy > ^max.Htf- ii) the 
intensity ratios between different events differ widely from line 
to line (e.g. H/5 is about two times more intense in #5 than in 
#6, Hy is only 0.8 times as intense in #5 as in #6. We have 
then also tried to fix the radial velocity of the gaussian H i com- 
ponents to the values obtained from the Fe ii lines in order to 
improve the fit, but the results were even worst: the method 
diverged for Hfi (i.e. no fit could be found) and the sequence 
of intensities for the Balmer series in event #6 was unphysical 

(^max,Hf = ^max.Hi = 2.5/?max,Hy, HO H/3 Component). 

One can also argue that the absence of Hi counterparts 
to the metallic events is a consequence of blending, which is 
unresolved by the multigaussian deconvolution procedure. In 
order to test such possibility, we have generated a composite 
synthetic R profile consisting of three gaussians with the same 
/^max and V as the Fell 5018 A JD 1023.52 spectrum and Av = 
150 km/s, which is the average velocity dispersion of the only 
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Fig. 13. Original spectrum (i.e. prior to the synthetic stellar 
spectrum subtraction) as observed in the JD 1023.52 spectrum 
of WW Vul. Wavelengths have been converted into radial ve- 
locities. Plotted lines are Hy, H6 and the three lines of the 
Fen 42 triplet. The photospheric synthetic spectrum used to 
form R is shown in dotted lines for Hy and Fe II 4924A (being 
very similar for the other lines). The synthetic spectra reveal 
the presence of emission in the wings of Hy and a little in the 
blue wing of Fe II 4924A. The average radial velocity <v> of 
the events identified after the photospheric spectrum removal 
and R multigaussian fit is indicated using the same colour and 
label conventions of Fig. |9l There is 1 blueshifted TAC (#6, 
blue colour) observed both in metallic and hydrogen lines, note 
that Ave I » Avpen. There are also 2 metalHc TACs (@5, @7, 
green) observed only in metallic lines. Event @7 can be clearly 
identified by visual inspection as a relatively narrow peak in 
the Fe n lines. In events #6 and @7 <v> does not exactly corre- 
spond with the local minima in the Fe ii spectrum because the 
whole R profile is used in the multigaussian fits. (This figure is 
available in color in electronic form). 



H I component detected in that date. We have also added ran- 
dom noise to achieve a S/N ratio worst than that of the Hfi and 
Hy lines of that spectrum; finally, we have applied our multi- 
gaussian deconvolution procedure to the synthetic profile. The 
result is that the three individual gaussians have been success- 
fully retrieved. Note that the gaussian parameters we have in- 
troduced are test values. This result suggests that we would 
have been able of identifying H i counterparts of the metallic 
kinematics components if they would exist. Therefore, we con- 
sider that the existence of metallic absorptions without obvious 
H I counterparts is a rather firm result. 

There are no obvious differences, in terms of kinematics 
and absorption strength, between the WW Vul metallic events 
and the rest of the events detected in this star and in BF Ori, 
SV Cep, XY Per and UX Ori. Some remarkable differences 
appear, however, when the ion column densities causing the 
events are compared. Lower limits on the ion column densities 
causing the absorptions can be estimated according to the fol- 



Table 5. Column density estimates (according to Eq.^ of ions 
excited to the energy ground level of B aimer and Fe ii 42 triplet 
Unes for events WW Vul #7 (metallic) and XY Per #8. 
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Hy 
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10'3 
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3.3 


10'2 
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10'3 
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2.9 


10'* 



lowing formula ('SD itzeill978l) . which becomes exact when the 
gas is optically thin: 

47reo nigC^WA 



(1) 



Where Ng is the column density of ions in the ground state 
of the line, e is the elementary charge, is the electron mass, 
eo is the permeability of free space, c is the speed of light, / 
is the oscillator strength of the line, A is the wavelength of the 
line and Wa is the equivalent width of the absorption compo- 
nent. The column densities obtained are similar, within an or- 
der of magnitude, for all the metallic events. An upper limit on 
the column density of the H i atoms in the B aimer energy level 
(n=2) can be estimated if we assume that TACs with an inten- 
sity lower than 3 times the noise level of the spectra cannot be 
detected. Table [S] gives the estimated values for the Feii lines 
and some Balmer lines in the case of event WW Vul #7. Event 
#8 in XY Per has extremely weak saturated Hi features and 
strong metallic components (H i line strengths are lower than 
or just comparable to the Feii lines); in this sense, this event is 
the more similar to the metallic ones among all detected events 
with Hi and metallic components. Table |5] gives lower limits 
to the column densities, estimated from the Balmer and Fe ii 
lines of event #8 in XY Per. The values of Table|5]nicely shows 
the difference between the estimated column densities of ions 
excited to the energy ground level of the Feii 42 triplet and 
Balmer lines in both type of events, and it could point out to a 
fundamental difference on the nature and origin of the gas from 
which the absorptions rise. However, a definitive statement on 
this issue requires a deep analysis of the metallic events de- 
tected in WW Vul in order to estimate chemical abundances, 
or abundance ratios among elements. Such analysis, which is 
beyond the scope of this paper, needs a NLTE treatment of 
the spectra; such treatment would produce model dependent re- 
sults, since a previous knowledge of physical quantities, as for 
example volume densities of the gaseous clumps and electron 
temperature, are required. Such quantities cannot be estimated 
from our data in a confident way. 

5. Concluding remarks 

We have analyzed optical high resolution spectra of the HAe 
stars BF Ori, SV Cep, WW Vul and XY Per These spectra 
monitor the stars on time scales of months and days and, as in 
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the case of the previously studied UX Ori (IPaper if), they pro- 
vide observational constraints, which should be considered for 
any realistic scenario of the gaseous circumstellar disks around 
intermediate-mass PMS stars. Our results and conclusions can 
be summarized as follows: 

1. The gaseous circumstellar environment of these stars is 
very complex and active. The spectra always show circum- 
stellar line absorptions with remarkable variations in their 
strength and dynamical properties. 

2. Variable absorption features are, in most cases, detected si- 
multaneously in hydrogen and in many metallic lines with 
similar velocities. In each case, there are several kinematic 
components in each line, both blue-shifted and red-shifted 
with respect to the systemic velocity, denoting the simulta- 
neous presence of infalling and outflowing gas. We attribute 
the variable features detected in both Balmer and metallic 
lines to gaseous clumps of solar-like composition, evolving 
dynamically in the circumstellar disks of these objects. In 
this respect, the disks around the stars studied in this paper 
ar e similar to the UX Ori disk. Following the conclusions 
of |Paper3 we suggest that these clumps and their dynam- 
ical evolution should be investigated in the context of de- 
tailed magnetospheric accretion models, similar to those of 
T Tauri stars. 

3. The star WW Vul is peculiar and behaves differently from 
the other stars studied in this paper and also from UX Ori. 
It is the only star that shows, in addition to events seen both 
in metallic and hydrogen lines, similar to those observed 
in the other stars, also transient absorption components in 
metallic lines that do not apparently have any obvious coun- 
terpart in the hydrogen lines. This result, taken at its face- 
value, would indicate the presence of a metal-rich gas com- 
ponent in the environment of WW Vul, possibly related to 
the evaporation of solid bodies. However, any such conclus- 
sion is premature. We think that a series of optical spectra 
with better time resolution (hours) and longer monitoring 
(up to around seven days), spectra in the far UV range - to 
analyze Lyman and metallic resonance lines - and detailed 
NLTE models of difi'erent CS gas environments are essen- 
tial for further progress and to provide clues on the origin 
of these apparently metal-rich events, in terms of their ap- 
pearance/disapperance statistics, dynamics, metallicity and 
nature. 
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Table 6. Identified Transient Absorption Components (TACs) in BF Ori. Column 1 gives the corresponding Balmer or metal- 
lic line, Col. 2 gives the Julian Date (-2450000), Col. 3 represents the event assigned to the particular absorptions (see text 
Sect. 13. 2L Cols. 4 to 6 give the parameters of each transient absorption as estimated from the multigaussian fit of the normalized 
residual absorption R: v, radial velocity, Av, FWHM, and /^max. the strength of the absorption (peak of the gaussian). "0" in Col. 
3 corresponds to the narrow IS absorptions, while "-" means that the absorption is not associated with a particular event. 
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Table 7. Identified Transient Absorption Components (TACs) in SV Cep. Details as for Table|6l 
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71 
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70 
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3 


68 
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77 
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3 


51 
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-3 


15 
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3 


70 
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1113.4730 





-1 


12 
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3 


60 
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-1 


16 
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4 


-6 
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-1 


14 


1.03 


Hr 
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2 


13 
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4 
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-1 


14 
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He 
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4 


15 
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_2 


23 
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H( 
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4 
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0.40 


H8 
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7 
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29 
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Hy 


1209.3372 


7 
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4 


-2 


19 
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4 


12 


33 
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He 


1209.3372 


7 
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298 


0.29 
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1025.6260 


4 


12 


39 


0.14 


H? 


1209.3372 


7 
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0.24 


Fen5169A 


1025.6260 


4 


10 


39 


0.15 


CanK 


1209.3372 


7 
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0.19 
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1025.6260 








17 


0.95 


Fen 4924A 


1209.3372 


7 
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0.05 
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-0 


14 


0.96 


Fen 5018A 


1209.3372 


7 
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0.07 


H/J 
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3 


80 
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0.75 
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1209.3372 


7 
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Hr 
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3 


72 


190 


0.76 
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8 


-4 


99 


0.97 
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3 


69 
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3 
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He 
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3 


87 
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0.56 
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3 


66 
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8 


6 
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3 


75 
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8 


9 


99 


0.46 


Nai D2 
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3 


82 
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0.16 














NaiDl 
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3 


75 
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0.12 
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1209.3372 


8 


40 


96 


0.07 


Fen 4924A 


1026.6684 


3 


84 


87 


0.11 


NaiDl 


1209.3372 


8 


33 


51 


0.06 
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1026.6684 


3 


87 


95 


0.16 
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1209.3372 


8 


19 


76 


0.09 
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3 


86 


83 
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8 


16 


78 
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H/J 


1026.6684 


4 


-8 


82 


0.59 
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1209.3372 


8 


17 


74 


0.15 


Hr 


1026.6684 


4 
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68 


0.42 


CanK 


1209.3372 





-I 


18 


0.43 


H6 
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4 
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60 


0.35 
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1209.3372 





1 


22 


0.77 


He 


1026.6684 


4 
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20 
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4 
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53 
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1 


15 
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Table 8. Identified Transient Absorption Components (TACs) in WW Vul. Details as for Table|6l 
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JD 
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1 
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95 


0.10 


Hr 


950.6176 


1 


71 
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0.41 


m 
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1 


58 
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0.44 


He 


950.6176 


1 


59 
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1 


69 
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1 


74 
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0.34 
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950.6176 


1 


40 
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1 


45 
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0.20 
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950.6176 


1 


44 


207 
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950.6176 


2 
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36 


0.17 
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950.6176 
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41 


0.22 
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2 
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41 
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3 
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-11 


58 
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NaiDl 


950.6176 


3 


-22 
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58 


0.63 


Hy 
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4 
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55 


0.43 
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4 
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55 


0.39 


H^ 
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4 
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49 


0.25 
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4 


-86 


42 


0.31 
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4 


-74 


55 
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4 


-67 


45 


0.09 
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4 
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46 


0.12 
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-4 
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1 
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0.33 
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1 


72 
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1 


42 
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0.68 
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1 


70 


91 


0.34 
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1 


30 
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0.57 
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1 


59 


62 


0.14 
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1 


67 


43 


0.12 


Fen 4924A 
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1 


45 
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0.18 
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1 


44 
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0.26 
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i 


46 
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951.6232 


2 
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53 
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951.6232 


2 
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46 


0.74 
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951.6232 


2 
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59 


0.50 
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951.6232 


2 


-10 


58 


0.36 
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951.6232 


2 
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42 


0.12 
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2 
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52 


0.21 
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951.6232 


2 
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66 


0.33 
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3 
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Hy 


951.6232 


3 
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3 


-53 
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He 
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3 
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0.79 
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3 
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52 


0.21 
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3 
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41 


0.31 
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3 
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50 
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3 
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3 
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3 
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84 
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4 
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63 


0.14 


Hy 
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4 
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81 


0.19 
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4 
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78 


0.18 


H^ 
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4 
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4 
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-3 
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28 
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50 
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8 


86 
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8 


96 
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1112.3689 


8 
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248 


0.11 


Tiu4444A 


1112.3689 


8 


91 
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0.03 


Tin4572A 


1112.3689 


8 


78 
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0.03 


H/J 


1112.3689 


9 


7 


156 


0.98 


Hy 


1112.3689 


9 


17 
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0.97 


H5 


1112.3689 


9 


17 
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0.99 


He 
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g 


24 
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93 


Can K 


1 1 12 3689 


9 


23 
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0.78 


Nai D2 


1 112.3689 


9 


21 


83 


0.59 


Nai Dl 


1 112.3689 


9 


22 


72 


0.45 


Fe n 4924A 
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9 


18 


77 
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17 
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9 


35 


79 


0.08 
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9 


17 


66 
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-7 


25 
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-5 


56 


0.89 
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-2 


14 
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-I 


13 


0.62 


Fe II 4924A 
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13 


0.11 
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1112.3689 





I 


12 
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1112.3689 





-0 


13 


0.15 


Sen 4247 A 


1112.3689 





2 


31 
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Hy 


1113.3958 


8 
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248 
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8 
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He 
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0.62 


*H 
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g 
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47 
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1 113.3958 


8 
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0.46 
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8 


82 
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8 
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8 


74 
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0.06 
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1113.3958 


8 
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0.06 
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1113.3958 


8 
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8 
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9 


6 


55 
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9 


6 


46 
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Table 9. Identified Transient Absorption Components (TACs) in XY Per Details as for Table|6l 
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